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We have studied interoccurrence time distributions by analyzing the synthetic and three natural 
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£\J ' catalogs of the Japan Meteorological Agency (JMA), the Southern California Earthquake Data 

Center (SCEDC), and Taiwan Central Weather Bureau (TCWB) and revealed the universal feature 
of the interoccurrence time statistics, Weibull - log Weibull transition. This transition reinforces the 
view that the interoccurrence time statistics possess Weibull statistics and log- Weibull statistics. 

■ 

' Here in this paper, the crossover magnitude from the superposition regime to the Weibull regime 

m\ is proportional to the plate velocity. In addition, we have found the region-independent relation, 
■ m 2 c /m max = 0.54 ± 0.004. 

O ' PACS numbers: 91.30.Dk, 91.30.Px, 05.65.+b, 05.45.Tp 
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C« 1 I. INTRODUCTION 
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1 Statistical properties of time intervals between successive earthquakes, henceforward the interoccurrence times and 
the recurrence times, have been frequently studied in order to predict when the next big earthquake will happen. 
Previous papers have been focused on the determination of the probability distribution and the presentation of the 



he determination ot the probability distribution and the pre 
Q, A S Q, Q, 3- For instance, the Weibull distribution Q], 



Q\ , scaling law, as shown in the works of [1|, [2], y, |4|, |5|, [6(. For instance, the Weibull distribution the exponential 
. distribution {2]] , the Brownian passage time (BPT) distribution jjj , the generalized gamma distribution ^, 0, Q| , and 
qq . the log normal distribution are candidates for the distribution function of interoccurrence times and recurrence 
■ times. Also, in the stationary regime, a unified scaling law was proposed by Corral Q and then obtained by analyzing 

the California aftershock data [6[ . 
. £h ' Meanwhile, the time interval statistics have been studied by numerical simulations of earthquake models Q, because 

real earthquake data are limited. For example, both the conceptual spring-block models l], ij] and the sophisticated 
?3 , Virtual California model [ljjl show the Weibull distribution of the recurrence times. One of the authors (TH) also 

reported that the survivor function of interoccurrence times in the 2D spring-block model can be described by the 

Zipf-Mandclbrot type power law ll|, which has been observed by Abe and Suzuki 12 1. 



Very recently, a statistical feature of the interoccurrence times, the Weibull - log Weibull transition, was proposed 
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TABLE I: List of Earthquake Catalogs 



Catalog Name 


Covering Region 


Term 


Number of Earthquakes 


Terrain 


m° c 


JMA [14] 


25°-50°N and 125°-150°E 


01/01/2001-31/10/2007 


130244 


2.0 


2.0 


SCEDC [17] 


32°-37°N and 114°-122°W 


01/01/2001-31/12/2007 


10838 


0.0 


2.0 


CWB [18] 


21°N-26°N and 119°-123°E 


01/01/2001-31/12/2007 


148155 


0.0 


2.0 


Synthetic 


50 x 50 (system size) 




297040 


0.0 


0.3 



by analyzing the Japan Meteorological Agency (JMA) catalog 13|. We found that the probability distribution of 
interoccurrence times can be very well fitted by the superposition of the log-Wcibull distribution and the Weibull 
distribution. In particular, the distribution of large earthquakes obeys the Weibull distribution with exponent less 
than unity indicating that the process of large earthquakes is long-term correlated. Our earlier results emphasize 
that the interoccurrence time statistics basically contain both Weibull and log- Weibull statistics, and as the threshold 
of magnitude m c is increased, the predominant distribution could change from the log- Weibull distribution to the 
Weibull distribution. Those statistical properties, including the Weibull - log Weibull transition, can be also found in 
synthetic catalogs produced by the 2D spring-block model ,gj]. However, the applicability to other tectonic regions of 
the Weibull - log Weibull transition remains unsolved. 

In this study, we further investigate the interoccurrence time statistics by analyzing the Southern California and 
Taiwan earthquake catalogs. Together with previous results from the JMA and synthetic catalogs, we show the 
universal Weibull - log Weibull transition obtained in all of the catalogs. We also demonstrate that a crossover 
magnitude, m 2 , between the superposition regime and the pure Weibull regime is proportional to the plate velocity, 
and at the end of this paper we elucidate its implication in the geophysical sense. 



II. DATA AND METHODOLOGY 



For studying the interoccurrence time statistics we analyzed the JMA 14 1 , the Southern California Earthquake Data 

18| . the synthetic catalogs. The information 



Center (SCEDC) [17J, the Taiwan Central Weather Bureau (TCWB) 
on each catalog is listed in Table [U m m i n and rrfi c correspond to the minimum magnitude and the minimum cutoff 
magnitude, respectively. We basically consider real earthquakes with magnitude greater than and equal to 2.0 because 
smaller earthquakes are strongly incomplete. 



The synthetic catalog is created by the 2D spring-block model with the velocity-weakening friction law 



15J. This 



model is characterized by five parameters: the stiffness l 2 and ly, the decrement of the friction force a, the plate velocity 
u, and the difference between the maximum friction force and dynamical friction force. We set those parameters at 
3.5,1/ = 0.01, and a = 0.01 for reproducing the realistic Gutenberg-Richter law ll|. Event 



I 2 = 1 I 2 



3, a 



magnitude m in the model is defined as m = mg + log 10 (V^j ^ u i,jj /l -5, where Su^j and n arc the total slip at 
and the number of slipping blocks, respectively, mo is set at 0.7, so as to shift m to positive. The occurrence time of 
an event is the time when a block slips for the first time during an event. 

The nth interoccurrence time r„ is defined as t n+ \ — t n , where t n and £ n +i are the occurrence times of the nth 
and n + 1th earthquakes, respectively. The methodology is the same as that used in our previous studies Q. 
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TABLE II: The interoccurrence time statistics on Southern California and Taiwan. Note that the rms value calculated from 
the cumulative distribution of P(r) is on order of [lCP 3 ]. (•) stands for the ratio of P w to P(t). 







Pw ~\~ P{ VJ 


Pw ~\~ Ppow 


Pyj ~\~ Pgam 


Pw ~\~ Pin 


Region 


m c 


rms (P w /P(t)) 


rms (P w /P(r)) 


rms (P w /P{r)) 


rms (P w /P(t)) 




2.0 


1.7 (0.45 ±0.006) 


7.1 (0.44 ±0.02) 


25 (0.97 ±0.05) 


15 (0) 




2.5 


2.3 (0.58 ±0.007) 


7.2 (0.72 ±0.01) 


17(1) 


6.0 (0) 


California 


3.0 


5.3 (0.79 ±0.02) 


6.8 (0.91 ±0.01) 


9.1 (1) 


5.4 (0.57 ±0.03) 




3.5 


7.6 (1) 


7.6 (1) 


7.6 (1) 


7.6 (1) 




4.0 


23 (1) 


23 (1) 


23 (0.34 ± 1.29) 


23 (1) 




2.5 


14 (0.40 ± 0.06) 


14 (0.20 ± 0.05) 


43 (0.96 ± 0.07) 


32 (0) 




3.0 


3.4 (0.49 ±0.02) 


7.5 (0.53 ±0.02) 


24 (0.94 ± 0.06) 


12 (0) 


Taiwan 


3.5 


5.0 (0.64 ±0.02) 


8.4 (0.74 ± 0.02) 


13 (0) 


6.1 (0) 




4.0 


3.5 (0.68 ±0.02) 


7.4 (0.89 ±0.01) 


9.3 (0.83 ±0.03) 


3.4 (0.40 ±0.03) 




4.5 


5.8 (1) 


5.8 (1) 


5.8 (1) 


5.8 (1) 




5.0 


12 (1) 


12 (1) 


12 (1) 


12 (1) 



Earthquakes with magnitude m above a threshold magnitude m c in a studied region were considered. Then, we 
calculated the interoccurrence times and regressed the interoccurrence times data in the time domain of r > h. It 
should be noted that the spatial division was carried out when we examined JMA data, h are set at 0.5 (days) and 
for real and synthetic earthquakes, respectively. 



Our previous works 



revealed that the probability distribution of interoccurrence times P(t) can be described 



by the superposition of the Wcibull distribution P w and the log-Wcibull distribution Pi w , namely, 



P(t)= P xP w + {1- P )xP 1v 



(1) 



where 



7 _i 
t \ 7 



■ cxp 



(log 7^) 



log(r/fc) 
logr^ 



7,r c ,<5, t' c , and k are constants, p is the ratio of P w divided by P(r). Obviously, P(t) is Pi w when p = and P w 
when p = 1. k is fixed at 0.5. The log-Wcibull distribution is constructed by the log modification of the cumulative 
distribution of the Weibull distribution. Unity S reduces the log-Weibull distribution to the power law distribution. 



III. RESULTS 



The interoccurrence time statistics with different m c for the SCEDC and TCWB data are studied. We use the 
Weibull distribution P w , log-Weibull distribution Pi Wl power law distribution P powl gamma distribution P ga m, log 
normal distribution Pi n and their superposition to fit the distribution of the calculated interoccurrence time data. As 
shown in Table. UH the distribution with larger m c fairly follows the Weibull distribution. Comparing the obtained 
rms value (please see explanation in the supplementary), we found that P(t) is best fitted by the superposition of the 
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log-Weibull regime Superposition regime Weibull regime 




50 100 150 200 0.5 1.0 1.5 2.0 2.5 3.0 3.5 10 20 30 40 50 60 70 

Interoccurrence time Interoccurrence time (day) Interoccurrence time (day) 



FIG. 1: The cumulative distribution of interoccurrence times for different catalogs, (a), (b), and (c) correspond to the log- 
Weibull regime, the superposition regime, and Weibull regime, respectively. For (b) and (c), the result of the synthetic data 
are demonstrated in the inset figures. 



TABLE III: The results of fitting parameters for different distribution functions. 



caption 


catalog 


m c 




7 




5 




P 


(a) 


Synthetic" 


0.7 


2.8 






5.64± 0.01 


51.3± 0.06 





(b) 


Synthetic" 


1.6 


2.8 


1.15± 0.005 


311± 0.73 


7.36± 0.04 


604± 2.04 


0.69±0.05 


(b) 


Japan 6 


2.5 


7.2 


2.91± 0.14 


0.79± 0.01 


1.10± 0.03 


1.43± 0.01 


0.31±0.02 


(b) 


California 


2.0 


5.7 


1.83± 0.02 


0.81± 0.006 


1.13± 0.01 


1.48± 0.006 


0.45±0.006 


(b) 


Taiwan 


2.5 


7.1 


3.35± 0.55 


0.60± 0.03 


1.01± 0.09 


1.18 ± 0.02 


0.40±0.06 


(c) 


Synthetic" 


1.6 


2.8 


1.29± 0.008 


1115± 3.05 






1 


(c) 


Japan 6 


4.5 


7.2 


0.82 ± 0.009 


7.69± 0.07 






1 


(c) 


California 


3.5 


5.7 


0.95± 0.005 


9.46 ± 0.04 






1 


(c) 


Taiwan 


4.5 


7.1 


0.92± 0.005 


5.44± 0.02 






1 



"Dynamical parameters are set at 1%. = 1, ly = 3, a = 3.5, v = 0.01, and a = 0.01. 

6 We analyzed the data in the region, ranging from 30°N to 35°N, and ranging from 140°E to 145°E. 



Weibull and log- Weibull distributions. The ratio of P w to P(t), p, increases as m c is increased. Thus, we can claim that 
a transition from the Weibull to the log-Weibull transition occurs in both earthquake catalogs of Southern California 
and Taiwan. As an example, together with our previous results for the JMA and synthetic catalogs [§, 13|, we show 
the cumulative distributions of the interoccurrence times exhibiting the log-Weibull regime, the superposition regime, 
and the Weibull regime in Figs.[T](a), (b), and (c), respectively. Note that Japan data (x) in Fig. Q] is obtained from 
a divided region spanning 35°-40°N and 140°-145°E. In addition, the pure log-Weibull regime can only be extracted 
from the synthetic data as shown in Fig. [1] (a). 

The maximum magnitudes, denoted by m max , for the fitting results of P(t) in Fig. [1] are listed in Table. Hill For 
each studied area 7 gradually decreases as m c is increased. In the Weibull regime, 7 obtained from real earthquakes 
is less than unity, whereas 7 for our synthetic catalog is greater than unity. However, we have noted that 7 derived 
from the synthetic data could be also less than 1 by tuning the parameters of the spring-block model. For all catalogs, 
as m c is increased, r c , S and t' c increases double exponentially, linearly and exponentially, respectively. Particularly, 
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FIG. 2: The relation between p and the scaled magnitude m^/m max . As for the synthetic data, the Weibull - log Weibull 
transition appears clearly. 



p gradually increases when m c is increased, indicating that the probability distribution changes as m c is varied. The 
interoccurrcncc time statistics exhibit the Weibull - log Weibull transition, which means that the dominant distribution 
changes from the log- Weibull distribution to the Weibull distribution with increase in m c . We thus define crossover 
magnitudes m\ and representing magnitudes that the distribution change from the pure log- Weibull regime to 
the superposition regime and from the superposition regime to the pure Weibull regime, respectively. We focus this 
paper on the values of w? c and they are 1.7, 4.0, 3.3, and 4.4 for the synthetic, Japan, California, and Taiwan catalogs, 
respectively. 

Most importantly, the universality of the Weibull - log Weibull transition can be clearly shown when we consider a 
scaled magnitude with m c /m max . Fig. [2] demonstrates the relation between p and m c /m maxi the scaled magnitude, 
for four catalogs. As shown in Fig. [2]p gradually increases with increases in rn c /m max , and the scaled crossover 
magnitude m^,/m max of the Weibull - log Weibull transition is estimated to be approximately 0.60 for all of four 
catalogs. 

To further discuss the feature of this transition, we summarize our results obtained from 17 different regions (14 
regions in Japan and other 3 regions of Southern California, Taiwan, and the synthetic data) in Fig.[3j We surprisingly 
find an area-independent constant for the scaled crossover magnitude, namely Trr c jm max = 0.54 ±0.004. Note that, in 
Fig. 3, there are three outliers whose values of m1/m max are less than 0.5. For those three regions of 30°-35°N and 
125°-130°E, and 45°-50°N and 140°-145°E, the total numbers of earthquakes are 1739 and 1135, respectively, which 
are smaller than the earthquake numbers used in other regions (see Fig. 1 in Ref. [13J). Therefore, the small values 
of rr? c jm max are probably caused by the statistical issue of insufficient samplings. As for the region of 30°-35°N and 
125°-130°E, m max is 8.0 which is the largest magnitude throughout the JMA catalog we analyzed. Therefore, the 
scaled crossover magnitude for this region tends to be biased small. 



6 



0.7 n 



0.6 



E 

E 0.5 
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o 25-30N and 125-1 30E 

o 25-30N and 130-1 35E 

a 25-30N and 140-1 45E 

x 30-35N and 125-1 30E 

3 + 30-35N and 130-1 35E 

o 30-35N and 130-1 35E 

□ 30-35N and 135-1 40E 
a 35-40N and 130-1 35E 
x 35-40N and 135-1 40E 
+ 35-40N and 140-1 45E 
o 40-45N and 135-1 40E 

o 40-45N and 140-1 45E x 

a 40-45N and 145-1 50E 

x 45-50N and 140-1 45E 

+ California 

o Taiwan 

□ Synthetic (BK model) 

_J I I 



Maxmum magnitude m max 

FIG. 3: The relation between the scaled crossover magnitude ml/m ma x and the maximum magnitude m ma x for different 
regions. ml ranges from 0.39 (45°-50°N and 140°-145°E) to 0.62 (Taiwan). 



TABLE IV: List of the cross-over magnitude and the plate velocity (see in Ref [20j, |2l|]). The notation of PH, EU, PA, and NA 
represent PHilippine sea plate, EUrasian plate, PAcific plate, and North American plate, respectively. 



Region 


relative plate motion 


velocity [mm/yr] 


ml 


Taiwan 


PH-EU 


71 


4.40 


East Japan 


PA-PH 


49 


3.80 a 


West Japan 


PH-EU 


47 


3.78 b 


California 


PA-NA 


47 


3.40 



"We take an average using three regions; 25°-30°N and 140°-145°E (mj = 3.7), 30°-35°N and 140°-145°E (m 2 c = 3.7), and 35°-40°N 
and 140°-145°E (mjl = 4.0). 

'We take an average using five regions; 25°-30°N and 125°-130°E (ml = 3.7), 25°-30°N and 130°-135°E (ml = 3.3), 30°-35°N and 
130°-135°E (m 2 c = 4.3), 30°-35°N and 135°-140°E (ml = 4.0), and 35°-40°N and 135°-140°E (ml = 3.6). 



IV. CONCLUSION AND DISCUSSION 



We have studied the interoccurrence time statistics of natural and synthetic earthquakes by analyzing the JMA, 
SCEDC, TCWB, and synthetic catalogs and found the universal Wcibull - log Weibull transition of the interoccur- 
rence distribution. We emphasize that interoccurrence time statistics contain both Weibull statistics and log- Weibull 
statistics. And, in this paper, we demonstrate the transition does occur for different tectonic settings. We also observe 
the area-independent scaling relation, namely rn^/m max = 0.54±0.004. Still the origin of the log- Weibull distribution 
and the Wcibull - log Wcibull transition remains open. Our present work represents the first step to fully understand 
the interoccurrence time statistics and the Wcibull - log- Wcibull transition for real earthquakes. 

Although the scaled crossover magnitudes m^./m max is area- independent, the crossover magnitude m 2 c from the 
superposition regime to the pure Weibull regime depends on the tectonic region. The geophysical implication for 
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the area-dependent m 2 c could be exposed when comparing the plate velocity with averaged m 2 .. As clearly seen 
from Table HV| m 2 c is on the average proportional to the plate velocity. That means the magnitude of the largest 
earthquake, m max , for a tectonic region is more or less proportional to the plate velocity since m^/0.54 = m max . Ruff 
and Kanamori showed a relation stating that the magnitude of characteristic earthquake which occurs in subduction- 
zone is proportional directly to the convergence rate 19(, which is thus consistent with our results. This study therefore 
suggests a possible physical interpretation for earlier observation about the velocity-dependence of the characteristic 
earthquake size. 
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